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ABSTRACT 


The  structure  of  the  flow  in  a  turbomachinery  passage  is 
extremely  complex.  Successful  turbomachine  design  requires 
detailed  information  about  the  flow  characteristics  in  order  to 
predict  the  performance  of  a  candidate  geometry.  The  designer 
needs  information  about  the  complex  interaction  between  the  blade 
rows  when  multistage  blade  rows  are  considered.  A  three- 
dimensional  steady  potential -flow  code  for  a  single  blade  row  has 
been  developed,  where  integral  equations  were  derived  to  model  the 
solid  surfaces  and  the  inlet  condition  for  both  external  and 
internal  flows  with  rotating  or  nonrotating  lifting  blades. 
Global  iterations  for  the  surface-panel  source  densities  and  the 
vortex  strengths  were  used  in  conjunction  with  a  Neumann  iteration 
for  solving  the  integral  equations.  In  this  paper,  two  numerical 
coupling  procedures  between  the  rotor  and  the  stator  are 
described.  The  phenomena  of  flow  interaction  between  the  rotor 
and  the  stator  are  studied  and  two  test  cases  are  presented.  The 
performance  of  a  vaneaxial  fan  is  predicted  using  this  procedure 
and  the  results  indicate  that  the  present  numerical  coupling 
procedure  can  be  used  by  designers  to  select  optimal  blade  shapes. 
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INTRODUCTION 

The  design  of  turbomachines  has  for  many  years  been  supported  by  an  ever 
increasing  volume  of  information  on  the  performance  aspects  of  the  machines.  This 
information  is  concerned  primarily  with  the  fluid  flow  problems  of  cascades  and  the 
meridional  flow  through  the  multiplicity  of  blade  rows  in  axial-flow  compressors  and 
turbines . 
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There  are  basically  two  recognized  forms  of  interaction  between  rotor  and  stator 
blade  rows.  One  is  due  to  the  wakes  behind  blades  or  other  elements  of  the  machine 
immersed  in  the  flow  due  to  viscous  effects,  and  the  other  is  due  to  potential -flow 
effects  which  would  occur  if  the  working  fluid  were  perfectly  inviscid.  All 
observations  made  in  practical  situations  represent  a  combination  of  the  two  sources 
which  it  has  not  been  possible  to  separate  except  on  the  basis  of  some  assumed 
knowledge  of  the  mechanisms  involved.  Reference  [1]  provides  an  example  of  an  axial 
flow  turbomachine  that  uses  both  upstream  inlet-guide-vanes  and  a  downstream  stator 
with  the  rotor.  The  blades  on  the  rotor  operated  at  high  relative  velocity  and 
could  therefore  be  expected  to  have  intense  wakes.  Removal  of  the  downstream  stator 
on  which  the  rotor  wakes  could  have  reacted  produced  a  slightly  measurable  reduction 
in  the  noise  level.  Removal  of  the  inlet -guide -vane ,  however,  produced  remarkable 
noise  reductions.  This  example  indicates  the  significance  of  the  potential -flow 
interaction  effect  between  the  inlet-guide-vane  and  the  rotor. 

Recently,  Navier-Stokes  solution  procedures  have  been  used  by  researchers  [2,3] 
in  seeking  viscous  solutions  for  interacting  flows  for  a  rotor -stator  operating  in 
the  compressible  range.  Detailed  flow  and  heat- transfer  information  can  be 
obtained.  Efforts  for  obtaining  these  solutions  are  tremendous.  Extension  of  these 
approaches  to  incompressible  flows  or  more  complex  geometries,  e.g.  asymmetric  blade 
row,  is  not  a  straight  forward  task  and  solution  strategies  have  not  been  practical. 

This  paper  uses  a  previously  developed  three-dimensional  steady  potential -flow 
code  [4]  for  a  single  rotating  or  nonrotating  blade  row  to  examine  the  potential- 
flow  effects  between  rotor  and  stator.  A  numerical  coupling  of  the  rotor- stator 
influence  is  developed.  A  mathematically  defined  centrifugal  machine  and  an  actual 
axial-flow  fan  were  used  to  demonstrate  the  present  scheme  for  predicting  machine 
performance. 
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THEORY 


Single -Blade -Row  Calculation 

A  potential -flow  solution  model  is  developed  to  investigate  the  rotor-stator 
interaction.  This  potential -flow  model  utilizes  a  surface -singularity  method  for 
the  representation  of  a  nonpenetrating  surface  condition,  and  an  inviscid  vortex¬ 
line  distribution  for  the  representation  of  a  lift- generation  blade  (4],  The 
surface  singularities  used  are  constant -strength  sources  or  sinks  on  each  panel. 
Thus  the  computed  flow  field  in  terms  of  the  strengths  of  the  surface  singularities 
(M)  and  the  vortices  (T)  is 

where  ^  is  a  scalar  potential  for  the  resultant  potential-flow  field,  M(Q)  and  T(Q) 
are  the  strengths  of  the  source  and  the  vortex,  and  rpQ  is  the  distance  between  a 
field  point  P  and  a  source/vortex  point  Q.  A  Fredholm  integral  equation  of  the 
second  kind  for  the  unknown  M(Q)  is  formulated  as 


2k  M(P)  -  f(P)  +  |  [  M(Q)  K(P,Q)  -  M(P)  K(Q,P)  ]  dS0 

+  [  2ir  -  o0(P)  ]  M(P),  (4) 

and 


8  1 

K(P,Q)  -  _  _ , 

3NP  rpQ 

<*o(P)  -  K(Q,P)  dSQ, 

JSo 

f(P)  -  v„(P)  -  V„(P), 


3 


where  K(Q,P)  is  the  transpose  of  the  kernel  of  K(P,Q),  N  is  the  surface  normal 
directed  into  the  flow  field,  V„  is  the  N-component  velocity  due  to  blade  rotation 
and  incoming  flow,  and  Vnr  is  the  N-component  velocity  due  to  vortices.  The  vortex 
strength  is  determined  by  using  the  Kutta  condition  near  the  trailing  edge  of  the 
blade.  In  this  paper  a  wake  tangency  condition  is  satisfied  at  the  Kutta  point. 

The  vortex  wake  structure  generally  adopts  a  helix  shape.  Its  modelling  for  moving 
and  nonmoving  blades  is  given  in  Appendix  A.  A  global  iteration  scheme  [4]  is 
developed  for  solving  both  the  unknown  strengths  of  sources  and  vortices. 

Rotor-Stator  Couplina 

Flow  past  an  axial -flow  or  a  centrifugal  single-stage  rotor  and  stator  is  truly 
unsteady  and  asymmetric.  The  circumferential  mean  of  the  flow  field  between  the 
blades  of  each  blade  row,  however,  can  be  represented  by  a  steady,  symmetric 
nonuniform  through  flow.  This  mean  flow  can  then  be  used  in  the  performance 
calculation.  The  single -blade -row  calculation,  described  in  the  previous  section, 
is  suitable  for  predicting  the  flow  field  with  this  symmetric  nonuniform  inlet 
condition. 

The  coupling  between  the  rotor  and  the  stator  can  be  modelled  by  considering 
either  the  stator  or  the  rotor  as  a  group  of  external  singularities  which  exert  an 
influence  on  the  other  one.  Therefore  when  the  flow  in  the  rotor  (or  stator)  is 
calculated,  it  experiences  an  effect  from  the  singularities  outside  the  rotor  (or 
stator).  To  include  this  effect,  the  first  term  on  the  right  hand  side  of  Eq.  (4) 
is  modified  as: 


f(P)  -  V„(P)  -  Vn,(P)  - 

VD*  (P), 


(5) 
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where  vm  is  the  N-component  velocity  due  to  external  singularities.  The  complete 
calculation  is  accomplished  by  an  iterative  scheme.  Initially  the  symmetric  flow  in 
the  rotor  is  calculated  using  the  single-blade-row  model,  then  the  asymmetric  stator 
flow,  due  to  relative  position  of  rotor  and  stator  passages,  is  calculated 
considering  the  calculated  rotor  singularities  as  a  set  of  external  singularities. 
After  completing  the  stator  flow  field  calculation,  the  first  iteration  is 
completed.  The  second  iteration  of  the  rotor  starts  with  the  influence  from  the 
external  singularities  of  the  stator  and  the  influence  from  the  now  asymmetric  flow 
between  the  rotor  and  the  stator.  A  similar  calculation  is  made  for  the  stator  flow 
field  to  complete  the  second  iteration.  This  iterative  procedure  can  be  carried  on 
until  a  converged  solution  is  obtained.  The  sources  and  vortex  of  each  blade  are 
adjusted  during  the  calculations  to  account  for  the  effect  of  the  asymmetric  flow. 
This  procedure  generally  requires  large  computer  memory  and  computing  time  for  each 
rotor  or  stator  calculation.  A  single  stage  centrifugal  pump  geometry  with 
mathematically  defined  blades  is  used  in  this  paper  to  demonstrate  this  approach. 

In  an  effort  to  reduce  the  computing  time  and  memory  required  for  the  previous 
approach,  a  technique  was  developed  that  takes  advantage  of  the  periodic  mean  flow 
field  through  the  rotor  and  the  stator.  After  the  initial  symmetric  rotor  solution 
is  obtained,  conservation  of  rotational  energy  is  preserved  between  the  rotor  exit 
flow  and  the  stator  inlet  flow.  A  preswirl  strength,  derived  in  Appendix  B,  is 
calculated.  This  preswirl  strength  is  added  to  a  symmetric  stator  flow  field  to 
preserve  the  energy  imparted  by  the  rotor.  Although  the  calculated  rotor  exit  flow 
is  nonuniform  in  the  tangential  direction  between  rotor  blades,  an  averaging  of  the 
calculated  velocities  is  performed  in  the  tangential  direction  at  each  radius. 
Therefore,  the  stator  inlet  flow  becomes  a  nonuniform  axial  flow  in  the  radial 
direction  combined  with  a  preswirl  strength.  These  inlet  conditions  enable  the 
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stator  flow  field  calculation  to  be  considered  as  axisymmetric .  For  the  second 
iteration,  the  external  singularity  effect  from  the  calculated  stator  solutions  is 
added  to  the  rotor  flow  field  calculations .  The  induced  velocity  on  each  rotor 
control  point  is  calculated  from  the  singularities  of  the  stator.  Then  the  boundary 
conditions  for  the  rotor  solutions  are  updated  by  Eq.(5).  Since  the  velocity 
contribution  from  the  external  singularities  in  Eq.(5)  decays  with  1/r2,  a  symmetric 
influence  from  the  stator  singularities  to  the  rotor  blade  row  is  assumed.  Results 
for  an  actual  one-stage  axial-flow  fan  using  this  simplified  procedure  are  given  in 
the  following  section. 
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RESULTS 


Mathematical  Blade  Rows 

In  Fig.  1,  the  geometry  of  the  impeller  and  the  diffuser  for  a  single  stage 
centrifugal  pump  is  defined.  The  five -blade  impeller  rotates  at  u  rpm,  but  the 
seven-blade  diffuser  is  stationary.  Simple  blade  shapes  with  zero  camber  for  both 
the  impeller  and  the  diffuser  are  selected  for  ease  of  geometry  generation  and 
numerical  calculation.  The  hub  and  shroud  are  horizontal  and  parallel  to  the  xy- 
plane.  There  is  no  inlet  pipe.  However,  the  inlet  swirling  effect  in  the  impeller 
entrance  is  included  by  the  prerotation  T0.  The  mesh  pattern  shown  in  Fig.  2  is  for 
one  complete  repeating  sector,  which  covers  60  degrees  for  the  impeller,  and  51.43 
degrees  for  the  diffuser,  on  the  hub.  There  are  938  and  688  panels  for  the  impeller 
and  diffuser  meshes,  respectively.  Due  to  the  resemblance  to  a  2-D  blade,  a  single- 
bound-vortex  model  without  trailing  vortices  is  used  here  for  simplicity.  The  bound 
vortices  are  located  at  x  -  0.425  ft  and  y  -  0.075  ft  for  the  impeller  and  x  -  0.7 
ft  and  y  -  0.225  ft  for  the  diffuser  in  Fig.  1.  The  Kutta  points  are  at  the  middle 
points  of  the  semicircles  which  form  the  blade  trailing  edges.  Since  an  asymmetric 
flow  calculation  is  adopted,  only  one  Kutta  point  is  used  for  each  impeller  or 
diffuser  blade.  The  calculated  relative  impeller  velocity  distributions,  for  the 
case  Vo*  -  0  in  Eq.(5),  are  shown  in  Fig.  3.  It  shows  the  distributions  for  low, 
near-design  and  high  flow  cases.  The  alignment  of  the  velocity  vectors  on  the 
suction  side  of  the  rotor  blade  indicates  that  flow  separation  will  occur  even  at 
the  "design  condition".  This  is  due  to  the  zero-camber  blade  used  and  would  not  be 
expected  in  a  design  with  properly  designed  blades.  At  low  flow  condition, 
pressure-side  separation  is  also  indicated. 

Three  iterations  were  performed  for  the  impeller-diffuser  interaction.  The  final 
computed  absolute  velocity  distribution  is  shown  in  Fig.  4.  The  convergence  of  the 
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solution  procedure  is  monitored  by  the  flow  continuity  calculations  at  both  the 
impeller  and  the  diffuser.  Figure  5  shows  the  fluid  mass  recovery  for  the  diffuser 
at  Q  -  1000  gpm.  Since  the  flow  control,  i.e.  inlet  boundary  condition,  is  at  the 
inlet  of  the  impeller,  the  flow  continuity  calculation  is  always  well  satisfied  for 
the  impeller.  The  continuity  calculations  for  the  diffuser  show  a  smooth  increase 
for  each  iteration  to  Q  -  1000  gpm.  The  diffuser  average  flows  at  three  different 
radii  for  each  iteration  are  about  76,  91,  and  97  percent  of  the  Impeller  prescribed 
inlet  flow.  The  variations  of  the  calculated  bound  vortices  and  flow  field  are 
similar  to  the  variations  of  the  diffuser  continuity  calculations  . 

The  fluid  torque  is  calculated  as  the  difference  between  the  fluid  angular 
momentum  at  the  outlet  and  the  inlet, 

Tf  -  p  Vrr0V8  dS  -  p  V^V,  dS,  (6) 

Jso  K 

where  Sn  and  S*  are  the  surfaces  at  the  outlet  and  the  inlet,  Vr  and  V#  are  the 
radial  and  tangential  velocity  components,  and  p  is  the  fluid  density.  The  head 
change  associated  with  the  fluid  torque  is  defined  as 

H  -  Tk  «  /  p  g  Q,  (7) 

where  g  is  the  gravitational  constant.  Figure  6  shows  the  characteristics 
calculated  according  to  Eqs.  (6)  and  (7).  The  circle  and  the  triangle  are  for  the 
symmetric  flow  through  the  impeller  at  I\>  -  0  and  -50  ft2/s«c*  The  cross  and  the 
square  are  results  from  the  impeller -diffuser  interaction  calculation  at  the 
impeller  speed  w  -  600  rpm,  Q  -  1000  gpm  and  T0  -  0.  The  cross  (x)  represents  only 
the  contribution  across  the  impeller.  The  square  is  for  the  combined  contribution 
of  both  the  impeller  and  the  diffuser.  The  results  show  that  the  fluid  torque 
through  the  impeller  for  the  impeller-diffuser  unit  is  the  same  as  that  for  the 
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Impeller  by  Itself.  Due  to  the  energy  transfer  from  the  kinetic  to  the  potential 
energy  in  the  diffuser,  the  overall  fluid  torque  for  the  unit  is  reduced. 

Axial -Flow  Fan 

A  ventilation  supply  fan,  with  a  17-blade  rotor  and  a  11-blade  stator  using  the 
computational  grid  shown  in  Fig.  7,  is  used  as  a  test  example  for  validating  the 
present  simplified  coupling  procedures.  Due  to  the  requirement  of  the  conservation 
of  the  rotational  energy  between  the  rotor  and  the  stator,  the  solution  converges 
rapidly.  The  continuity  calculation  indicates  two  iterations  are  adequate  when 
overall  flow  quantities  are  concerned.  After  the  second  iteration,  the  rotor  and 
the  stator  velocity  distributions  at  u  -  3550  rpm  and  Q  -  16,000  cfm  are  shown  in 
Figs.  8  and  9,  where  the  plotting  locations  for  the  velocity  distributions  are  close 
to  the  blade  trailing  edges.  These  calculations  were  performed  with  the  trailing 
vortex  model  described  in  Appendix  A.  The  calculated  vortex  strengths  across  the 
spans  of  the  rotor  or  the  stator  blades  are  very  uniform.  Due  to  the  cancelation 
between  the  upper  and  lower  horseshoe  vortices,  the  results  indicate  that  the  blade 
flow  at  midspan  can  be  closely  modeled  by  a  two-dimensional  simulation.  The  tip  and 
hub  vortices  are  not  cancelled  and  have  a  dominant  effect  on  the  calculated  results. 
Effects  from  the  tip  and  root  vortices  are  shown  in  Figs.  8  and  9.  Figure  10  shows 
the  calculated  pressure  distributions  on  the  rotor  and  stator  blades.  Again  the  tip 
and  root  vortices  are  responsible  for  the  local  jump  in  the  pressure  distributions. 
For  other  parts  of  the  blades,  the  pressure  distributions  are  generally  uniform. 

The  overall  efficiency,  defined  as  the  ratio  of  the  rotor  energy  transfer  rate  to 
the  input  electric  energy  (IHP),  is  calculated  from  the  performance  prediction  and 
compared  to  the  measured  overall  efficiency.  Figure  11  shows  a  comparison  between 
the  calculated  efficiency  and  the  measured  performance,  which  has  a  measurement 
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uncertainty  of  ±1.2%.  Dependency  of  the  numerical  solution  with  grid  size  was  given 
in  reference  [4]  for  the  present  model.  Although  the  local  calculated  blade  surface 
pressure  between  fine  and  coarse  grids  can  vary  up  to  10%,  the  overall  efficiency 
from  the  numerical  results  is  estimated  less  than  4%.  The  comparison  shown  in  Fig. 
11  indicates  that  the  loss  due  to  the  fluid  energy  transfer  is  small. 


CONCLUSIONS 

A  computational  scheme  for  coupling  the  flow  field  predicted  by  a  potential -flow 
solution  method  for  a  rotor-stator  combination  has  been  developed.  Two  different 
coupling  procedures  have  been  evaluated:  asymmetric  and  symmetric.  The  convergence 
of  the  asymmetric  coupling  approach  was  demonstrated  by  an  impeller -diffuser 
geometry  with  mathematically  defined  blades.  The  symmetric  averaging  coupling 
approach  was  verified  using  an  actual  axial -flow  fan.  The  predicted  overall  machine 
performance  agrees  well  with  the  measured  performance  of  the  axial -flow  fan.  In 
conclusion,  it  has  been  shown  that  the  simplified  averaging  approach  can  be  used 
when  overall  averaged  quantities  of  the  flow  field  are  sought.  This  technique  is  an 
important  step  in  the  development  of  the  capability  to  design  and  analyze 
turbomachinery  designs  without  having  to  fabricate  and  test  each  iteration.  The 
result  of  this  will  be  to  increase  the  success  rate  while  at  the  same  time  reducing 
the  time  and  the  cost  to  develop  advanced  turbomachinery  for  the  Fleet. 
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APPENDIX  A 


Modified  Wake  Vortex  Model 

The  helix- shape  trailing  vortex  model  used  in  [4]  for  an  axial -flow  machine  is 
given  by  Eq.(15)  of  [4], 

x  -  rt  sin  8,  y  -  rt  cos  9,  z  -  zt  +  V.  8  /  u,  (Al) 

where  the  subscript  t  refers  to  the  trailing  edge,  and  8  is  the  azimuthal  angle 
measured  from  the  trailing  edge.  Equation  (Al)  can  also  be  written  as 

8  -  Az  w  /  V.,  (A2) 

where  Az  is  the  axial  distance  measured  from  t.  This  equation  indicates  that  the 
helix  is  independent  of  the  blade  shape  or  inlet  flow  angle.  A  modified  wake  vortex 
model,  which  is  suitable  for  moving  as  well  as  nonmoving  blades,  is  used  in  the 
present  calculation ,  i . e . 

rt  8  /  Az  -  (  rt  w  +  V«^t  )  /  V^t,  (A3) 

where  subscripts  A  and  T  represent  axial  and  tangential  directions.  Physically, 
this  model  relates  the  trailing  vortex  to  the  blade  trailing  edge  velocities.  There 
are  two  options  added  to  this  modified  relation.  First,  a  fixed  flow  exit  angle  at 
t  can  be  specified.  A  relaxation  parameter  can  be  input  to  control  the  vortex 
recovery  distance  back  to  the  state  of  incoming  flow. 
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APPENDIX  B 


Stator  Preswirl  Model 

The  prerotation  effect  can  be  simulated  by  a  line  vortex  along  the  axial 
direction  from  z  -  0  to  z  -  <*>.  If  the  strength  of  this  vortex  is  assumed  to  be 
constant  r0>  the  tangential  velocity  at  P  due  to  Te  is 


V,  - 


4  it  rp 


[  1  + 


Zp 


(  rPz  +  zPz  ) 


rr r 


]. 


(Bl) 


where  rP  is  the  radius  of  point  P  from  the  center  axis.  When  the  conservation  of 
rotational  energy  is  applied  to  the  rotor  exit  plane  and  the  stator  inlet  plane, 


JVX  rP  V#  dSP  -  V,  rP  V9  dSP. 

So  K 

Substituting  Eq.(Bl)  into  Eq.(B2),  one  obtains 


(B2) 


T0  -  4ir 


I, 


Vz  rP  Vg  dSP 


(B3) 
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30  ft/sec 


Figure  3.  Relative  Velocity  Distributions  at  Midspan  for  the  Synmetric 


Flow  Past  the  lapel ler,  «  -  600  rpm,  r.  -  0. 

(A)  Q  -  400  gpa,  (B)  Q  -  2500  gpa,  (C)  Q  -  4400  gpa. 


30  ft/sec 


Figure  4.  Absolute  Velocity  Distribution  at  Midspan  for  Flow  Past  the 
Impeller  and  the  Diffuser  After  the  Third  Iteration, 
u  -  600  rpm,  Q  -  1000  gpm,  r,  -  0. 


Flow  (GPM  x  10 


Iteration  No 


Figure  5.  Flow  Continuity  Calculations  at  Different  Radii  Between 
Successive  Iterations  for  Q  -  1000  gpm. 
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Figure  6.  The  Characteristics  for  the  Iopeller-Dif fuser 
(A)  Torque,  (B)  Head. 
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Figure  7 .  Computational  Grid  for  the  Air  Supply  Fan 
(Outside  casing  is  not  shown) 
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Figure  9.  Stator  Velocity  Distributions  near  Trailing  Edge, 

-  3550  rpm,  Q  -  16,000  cfm,  (A)  Axial  (B)  Radial  and 
Tangential  Components. 
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Figure  10.  Calculated  Blade  Pressure  Distributions  for  Single-Stage 
Ventilation  Supply  Fan,  u  -  3550  rpn,  Q  «■  16,000  cfm, 

(A)  Rotor  (B)  Stator. 
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Figure  11.  Performance  Prediction  for  Single-Stage* Ventilation  Supply 
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